Numerical integration of ordinary differential equations on manifolds

Summary This paper is concerned with the problem of developing numerical integration algorithms for differential equations that, when viewed as equations in some Euclidean space, naturally evolve on some embedded submanifold. It is desired to construct algorithms whose iterates also evolve on the same manifold. These algorithms can therefore be viewed as integrating ordinary differential equations on manifolds. The basic method “decouples” the computation of flows on the submanifold from the numerical integration process. It is shown that two classes of single-step and multistep algorithms can be posed and analyzed theoretically, using the concept of “freezing” the coefficients of differential operators obtained from the defining vector field. Explicit third-order algorithms are derived, with additional equations augmenting those of their classical counterparts, obtained from “obstructions” defined by nonvanishing Lie brackets.     

On the Moments of Overflow and Freed Carried Traffic for the GI/M/C/0 System

In circuit-switched networks call streams are characterized by their mean and peakedness (two-moment method). The GI/M/C/0 system is used to model a single link, where the GI-stream is determined by fitting moments appropriately. For the moments of the overflow traffic of a GI/M/C/0 system there are efficient numerical algorithms available. However, for the moments of the freed carried traffic, defined as the moments of a virtual link of infinite capacity to which the process of calls accepted by the link (carried arrival process) is virtually directed and where the virtual calls get fresh exponential i.i.d. holding times, only complex numerical algorithms are available. This is the reason why the concept of the freed carried traffic is not used. The main result of this paper is a numerically stable and efficient algorithm for computing the moments of freed carried traffic, in particular an explicit formula for its peakedness. This result offers a unified handling of both overflow and carried traffics in networks. Furthermore, some refined characteristics for the overflow and freed carried streams are derived.     

利用杨辉三角形对称性推导高阶运动微分方程

利用Mathematica数学软件计算函数r=r(q(t),t)各变量之间偏导和高阶导数的关系，发现具有杨辉三角形对称性.结合杨辉三角形的对称性规律和牛顿第二定律推导出了高阶运动微分方程，并讨论了理想约束系统下的高阶运动微分方程. 关键词： 杨辉三角形 牛顿第二定律 高阶运动微分方程 高阶力变率 高阶速度能量 理想约束     

SINGULAR PERTURBATION OF BOUNDARY VALUE PROBLEM FOR NONLINEAR HIGHER ORDER ORDINARY DIFFERENTIAL EQUATIONS INVOLVING TWO SMALL PARAMETERS

The singularly perturbed boundary value problem for nonlinear higher order ordinary differential equation involving two small parameters has been considered. Under appropriate assumptions, for the three cases:ε/μ2→0(μ→0),μ2/ε→0 (ε→0) andε=μ2, the uniformly valid asymptotic solution is obtained by using the expansion method of two small parameters and the theory of differential inequality.     

Attractors for a Three-Dimensional Thermo-Mechanical Model of Shape Memory Alloys

Abstract In this note, we consider a Frémond model of shape memory alloys. Let us imagine a piece of a shape memory alloy which is fixed on one part of its boundary, and assume that forcing terms, e.g., heat sources and external stress on the remaining part of its boundary, converge to some time-independent functions, in appropriate senses, as time goes to infinity. Under the above assumption, we shall discuss the asymptotic stability for the dynamical system from the viewpoint of the global attractor. More precisely, we generalize the paper [12] dealing with the one-dimensional case. First, we show the existence of the global attractor for the limiting autonomous dynamical system; then we characterize the asymptotic stability for the non-autonomous case by the limiting global attractor. * Project supported by the MIUR-COFIN 2004 research program on “Mathematical Modelling and Analysis of Free Boundary Problems”.     

The Neumann problem for some degenerate elliptic equations

In the paper we study the equation L _u = f, where L is a degenerate elliptic operator, with Neumann boundary condition in a bounded open set μ. We prove existence and uniqueness of solutions in the space H(μ) for the Neumann problem.     

EXISTENCE OF PERIODIC SOLUTIONS OF SYSTEM OF DIFFERENCE EQUATIONS

This paper studies the nonautonomous nonlinear system of difference equationsΔx(n)=A(n)x(n)+f(n,x(n)),n∈Z,(*) where x(n)∈R~N,A(n)=(a_(ij)(n))N×N is an N×N matrix,with a-(ij)∈C(R,R) for i,j= 1,2,3,...,N,and f=(f_1,f_2,...,f_N)~T∈C(R×R~N,R~N),satisfying A(t+ω)=A(t),f(t+ω,z)=f(t,z) for any t∈R,(t,z)∈R×R~N andωis a positive integer.Sufficient conditions for the existence ofω-periodic solutions to equations (*) are obtained.     

共振下n维Lienard型方程的2π周期解

本文在跨共振点条件下证明n维Lienard型方程存在2π周期解.     

A possible counterexample to well posedness of entropy solutions and to Godunov scheme convergence

A particular case of initial data for the two-dimensional Euler equations is studied numerically. The results show that the Godunov method does not always converge to the physical solution, at least not on feasible grids. Moreover, they suggest that entropy solutions (in the weak entropy inequality sense) are not well posed.

     

冻干溶液的低温显微研究与热分析

利用低温显微系统研究了叔丁醇水溶液降温时的晶体形态,并利用差示扫描量热法(DSC)研究了溶液的反玻璃化结晶。低温显微实验表明:叔丁醇／水溶液结晶形状、大小与叔丁醇的浓度有关,浓度为10％时形成粗大的针状结晶, 有利于升华速率的提高。DSC实验表明:10％叔丁醇／10％蔗糖／水三元溶液的临界复温速率为250℃／min左右。